Tight and loose shapes in flat entangled dense polymers 
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We investigate the effects of topological constraints (entanglements) on two dimensional polymer 
loops in the dense phase, and at the collapse transition (O point). Previous studies have shown that 
in the dilute phase the entangled region becomes tight, and is thus localised on a small portion of 
the polymer. We find that the entropic force favouring tightness is considerably weaker in dense 
polymers. While the simple figure-eight structure, created by a single crossing in the polymer loop, 
localises weakly, the trefoil knot and all other prime knots are loosely spread out over the entire 
chain. In both the dense and conditions, the uncontracted knot configuration is the most likely 
shape within a scaling analysis. By contrast, a strongly localised figure-eight is the most likely shape 
for dilute prime knots. Our findings are compared to recent simulations. 

PACS numbers: 87.15.-v, 82.35.-x, 02.10.Kn 



I. INTRODUCTION 



Knots, and topological constraints in general, play an 
important role in macromolecular systems. In gels and 
rubbers, permanent entanglements strongly influence the 
equilibrium and relaxation properties 0^ 0, 0, 01 ■ Even 
single molecules with identical chemical structure but dif- 
ferent topology may exhibit different physical properties 
Knots are also present in biological molecules: For 
example, some proteins exhibit knotted configurations 
0, and the active modification of DNA knots through 
energy-consuming enzymes (topoisomerases) poses inter- 
esting challenges to the issue of the knot detection 0, ■ 
Experimentally, the observation of individual molecules 
by single molecule force spectroscopy has come of age 
[3, ; these methods can be used to probe the mechan- 
ical behaviour of knotted biopolymers directly. 

Given a knot in a closed ring, an obvious question is 
whether, on average, the knot segregates into a small re- 
gion in which all topological details are confined, and a 
large, simply connected segment; or whether it is loosely 
spread over the entire chain. In the dilute phase, it has 
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been found by numerical evidence that flattened (hence 
two-dimensional [Tll | ) k nots are localised, i.e. tight [T^ . 
In a previous study we developed an analytical ap- 
proach, based on scaling results for polymer networks 
[14 [13, which explains and quantifies the tightness of 
any prime knot in fiat (2D) dilute polymer loops. Some 
aspects of our results have been verified experimentally 
by means of a vibrated granular chain It should be 
mentioned that DNA chains have experimentally been 
flattened by adsorption on an adhesive membrane |l7j . 
In 3D, the global topological constraints of a knot are 
hard to implement by analytical methods sim- 
pler topological invariant, the linking number of closed 
DNA rings, can be incorporated in the mapping to a 
field theory 0| by means of suitable gauge fields (l3|.) 
Consequently, the tightness or localisation of 3D knots 
has not been conclusively characterised. A number of 
phenomenological models and numerical studies support 
the localisation of simple knots |23| • Tightness has also 
been found in 3D slip-linked polymer chains in the dilute 
phase 21] . Conversely, delocalisation has been predicted 
for more complicated knots [2^ . 

In many situations, however, polymer chains are not 
dilute. Polymer melts, gels, or rubbers exhibit fairly high 
densities of chains, and the behaviour of an individual 
chain in such systems is significantly different compared 
to the dilute phase 0, Q . Similar considerations apply to 
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biomolecules: In bacteria, the gyration radius of the al- 
most freely floating ring DNA is often larger than the cell 
radius itself. Moreover, under certain conditions, there is 
a non-negligible osmotic pressure due to vicinal layers of 
protein molecules, which tends to confine the DNA [23l| . 
In protein folding studies, globular proteins in their na- 
tive state are often modelled as compact polymers on a 
lattice (see plj for a recent review). 

Given this motivation, in this work we consider self- 
avoiding dense polymers with permanent entanglements, 
complementing our previous studies for the dilute phase 
[i3,|21|. A polymer is considered dense if, on a lattice, 
the fraction / of occupied sites has a finite value / > 0. 
This can be realised by a single polymer of total length 
L inside a box of volume V and taking the limit L oo, 
y ^ oo in such a way that / = L/V stays finite |2^ 
|2^ . Alternatively, dense polymers can be obtained in an 
infinite volume through the action of an attractive force 
between monomers. Then, for temperatures T below the 
transition (Theta) temperature 0, the polymers collapse 
to a dense phase, with a density / > 0, which is a function 
of T (26> .27ij . For a dense polymer in d dimensions, the 
exponent defined by the radius of gyration Rg ^ L'^ , 
is simply v — 1/d. In 2D, the dense polymer phase for 
< / < 1 is related to a conformal field the ory, and the 
resulting scaling behaviour is known exactly |2^ |2^ . 
The limit / = 1 is realised in Hamiltonian paths, where 
a random walk visits every site of a given lattice exactly 
once. For some ca ses, their scaling behaviour is known 
exactly as well [23, 1^3 • Another way to make polymers 
collapse in 2D is to exert a pressure on self-avoiding loops 
(2D vesicles), conjugate to their area, which results in 
double- walled, branched structures _31, 32J. Recently, 
the corresponding crossover scaling function from linear 
to branched polymers (lattice animals) has been obtained 
exactly [s^. It is believed that the branched and dense 
polymer phases are different. 

Here, we extend our previous scaling analysis of knots 
in dilute polymers 13] to the dense phase, and at the 
Q point. The general conclusion is that the entropic 
mechanism for tightening of entanglements is consider- 
ably weaker in the denser regimes. While the simple 
figure-eight structure is still tight in these phases, the 
trefoil becomes loose, with the trends more pronounced 
in the dense phase. Note that at the O point in 2D, the 
swelling exponent is v = 4/7, implying a fractal dimen- 
sion l/v smaller than d = 2, and an asymptotic density 
of / = 0. In a recent numerical study, Orlandini, Stella, 
and Vanderzande (s^ (hereafter referred to as OSV) in- 
vestigate the tightness of the 2D projection of the trefoil 
knot, and find delocalisation, in contrast to the strong 
localisation obtained in the dilute case ^3 ■ 

The main focus of this work is thus the localisation of 
flat entangled polymers in the dense phase and at the 
@ point by means of scaling arguments, in analogy to 
the dilute case 0|. When possible, we compare these 
results with our own numerical simulations, as well as 
those by OSV. In section ^1 we first review the differen- 



tiation between tight (localised) and loose (delocalised) 
segments in entangled polymers. In section UTTl we then 
consider the 2D figure-eight structure (F8) and compare 
our scaling results with Monte Carlo simulations. In sec- 
tion llVl we derive the scaling results for the 2D projection 
of the trefoil, and compare them with the simulations by 
OSV. Section Ivl contains our conclusions. In appendix 
IXI we compile some scaling results for general polymer 
networks. Finally, at the end of appendix ^ we consider 
briefly a newly discovered phase of 2D dense polymers 
where the strict non-crossing condition is relaxed |28|. 



II. TIGHT AND LOOSE SEGMENTS 

Consider an entangled polymer chain in 2D of length L, 
such as a simple flat, once- twisted ring with one crossing 
(called "figure-eight" (F8), see position III in figure 
top row) or the 2D projection of the trefoil pressed flat 
against a surface by an external force (figure The 
orientations of the crossings are irrelevant, and can thus 
be considered as vertices with 4 outgoing legs. Thus, 
each structure is mapped on a 2D polymer network with 
a number of vertices which are joined by Af segments of 
variable lengths {si} under the constraint X]i=i ~ ^■ 
In the following we shall use the convention si < S2 < 

... <S^f. 

Since we are interested in the tightness of such a 
network, we define the size of the entangled region as 
£ — X^i^T^ so t^^t remaining (largest) segment is 
s^r = L - £. Clearly, for the F8, A/" = 2 and i ^ si, 
while for the trefoil J\f = 6 and £ = ^i- Note that 

the above definition does not necessarily imply that £ is 
small; however, if £ ^ L, the structure assumes the form 
of a possibly multiply connected knot region of size £ and 
a large simple loop of size L — £. For a knot represented 
by the network Q, an important quantity is its number 
of configurations ujg{£,L) for fixed £. In general, in the 
tight limit i L, the configuration number scales as a 
power law 



Log{£,L) ^ £-^^ [£-€.L) 



(1) 



where the exponent cg depends on the topology of the 
network Q 37]. 

For given lo{£, L), various quantities of interest can be 
calculated. For example, the mean size of the knot region 
is given by 



f-L/2 

{£) = / d££p{£,L) 



(2) 



with the (normalised) probability density function (PDF) 



p{£,L) = Loi£,L) 



L/2 



d£uj{£, L) 



(3) 



We have introduced a short-distance cutoff a, set by the 
lattice constant, to control the non-integrable singularity 



Dense 




0.125 



0.75 



1.375 




1.875 



6-point 






0.714 



1.143 



1.571 



2.095 



2.524 3.857 



Dilute 







2.6875 



3.375 



4.0625 5.5 



6.1875 



9.0625 



Figure 1: "Contractions" of tlie flat trefoil knot, arranged according to higher scaling orders. The protruding legs stand for the 
remaining simply connected ring. Below the contractions, we list the associated scaling exponent c (see text). Both for dense 
polymers and polymers at the point, the leading shape is the original (uncontracted) trefoil configuration on the very left. 
For comparison, we also show the corresponding ordering for the dilute phase 

at ^ = which occurs in equation ^ for c > 2, and in the denominator of equation (j^J for c > 1. 
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Thus, depending on the value of c, three cases can be 
distinguished: 

(i) c < 1: Both integrals are well-defined for a ~ 0, 
and one immediately obtains {£) ~ L due to dimensional 
considerations. The knot region grows linearly with L 
and is thus delocalised, i.e. spanning the whole polymer. 

(ii) 1 < c < 2: To leading order in a, one finds {£) ^ 
a'^~^L^^^. The size of the knot region scales with L, but 
with an exponent 2 — c < 1, and is thus weakly localised, 
i.e., {tj/L -> for L ^ oo. 

(iii) c > 2: To leading order in a, one finds {tj ^ a, 
i.e., the size of the knot region is independent of L for 
L oo, and is thus strongly localised. 

In a more complicated network, such as for the tre- 
foil in figure ^ several segments can be simultaneously 
tight. We shall refer to the emerging structures as possi- 
ble shapes (or contractions) of the original network. Each 
shape corresponds again to a network Q. The parameter 
c not only determines the tightness of the knot region for 
given shape Q, but also controls the overall likelihood Pg, 



of different possible shapes. The latter is obtained (up 
to normalisation) by integrating u!g{£,L) over all possi- 
ble values of £, as in the denominator of equation 
Depending on the values of c, different cases can be dis- 
tinguished: 

(i) c < 1: The integral is convergent for small £ with a 
finite limit for a — > 0, and scales as {a/Ly~^^. If there is 
a variety of shapes with c < 1, the one with the lowest 
value Cm is the most likely, Pg for the others scaling with 
a factor of (a/L)'^"'^'". This is expected, since networks 
with c < 1 are delocalised over the whole chain of length 
L. 

(ii) c > 1: The integral is divergent for small £, and is 
thus dominated by the lower cutoff. To make sure that 
the knot region is large enough to consider different ver- 
tices (crossings) as separated on a scale much larger than 
a, one should use a lower cutoff A of intermediate length 
a ^ A <^ L [s^- The relative probabilities for different 
shapes with c > 1 now scale with a factor of {a/ Ay. Dif- 
ferent shapes with c > 1 scale differently with the lower 
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cutoff A, but not differently witfi L. Tliis reflects tlie 
fact that for c > 1 tlie probability Pg is dominated by 
the small i behaviour of LOg{£, L). 

In both cases Pg scales with a factor of a"^. This is 
the scaling order of such a shape Q, as the likelihood of 
possible shapes with distinct values of c can be ordered 
according to powers of a'^ with more dominant networks 
for smaller values of c [l^ . 

III. THE 2D FIGURE-EIGHT: SCALING 
ANALYSIS & SIMULATIONS 

The most elementary entangled object in 2D is the F8, 
which consists of two loops of variable lengths £ and L — £ 
(position III in figure ^ top row). The crossing point can 
be considered as a vertex with 4 outgoing legs, and the 
F8 corresponds to a network with M = 2 segments of 
lengths si — i and S2 = L — i The number of 

configurations u!s{£,L) for the F8 can be deduced from 
results for general polymer networks, obtained by Du- 
plantier and coworkers 25, 26, 29,. 35J, which are com- 
piled in appendix IXI In 2D it has the scaling form 

L) ~ iv,{L) {L - ey^X (^-^) . (4) 

As discussed in the previous section, the localisation of 
F8 is controlled by the limit of ujg{£,L) for £ ^ 0, i.e., 
the behaviour of the scaling function X{x) for x — £/{L — 
£) 0. The latter can be determined by the following 
argument (see text below equation (59) in reference [2^'): 
Clearly, for £ <§; L, the big loop of length L — t will behave 
like a simple ring, so that los{£, L) should reduce to ujq{L). 
This implies Xix) ~ x'^'^ as x ^ and thus 

u;^{£,L)^u;,,{L)£-- {£<^L). (5) 

The value of the exponent c for the 2D dense F8 (see 
appendix is 

c = -78 = 11/8 = 1.375 , (6) 

implying that the smaller loop is weakly localised. This 
means that the probability for the size of each loop is 
peaked at € = and, by symmetry, at £ = L. An analo- 
gous reasoning for the 2D F8 at the point gives 

c= 11/7= 1.571 , (7) 

i.e., in this case the smaller loop is also weakly localised. 

Figure|2]shows the symmetric initial and a typical equi- 
librium configuration for periodic boundary conditions 
obtained from Monte Carlo (MC) simulations, see be- 
low for details. In figure 13 the lines represent the bonds 
(tethers) between the monomers (beads, not shown here). 
The three black dots mark the locations of the tethered 
beads forming the slip-link in 2D, by which we model 
the crossing jl3. |2]| . The initial symmetric configuration 



soon gives way to a configuration with £ <C i on ap- 
proaching equilibrium. Figure |21 shows the development 
of this symmetry breaking as a function of the number 
of MC steps. We note, however, that the fiuctuations of 
the loop sizes in the "stationary" regime appear to be 
larger in comparison to the dilute case studied in refer- 
ence |23|. We checked that for densities (area coverage) 
above 40% the scaling behaviour becomes independent of 
the density. The above simulation results correspond to 
a density of 55%, which is roughly half of the maximal 
possible density of 90% (closed packed area coverage). 




Figure 2: Symmetric {£ — L/2 = 128) initial configuration of 
a 2D dense F8 (left), and its equilibrium configuration (right) 
with periodic boundary conditions. The slip-link is repre- 
sented by the three (tethered) black dots. 



600 




Monte Carlo steps 

Figure 3: The breaking of the initial symmetry between the 
two loops of the F8, as a function of MC steps at 55% area 
coverage. 

The size distribution data is well fitted to a power law 
(for over 1.5 decades with 1024 monomers), and the cor- 
responding exponent with 512 and 1024 monomers in fig- 
ure 0] is in good agreement with the predicted value. 

For our MC analysis, we used a hard core bead-and- 
tether chain, in which self-crossings were prevented by 
keeping a maximum bead-to-bead distance of 1.9 times 
the bead diameter, and a maximum step length of 0.27 
times the bead diameter (cf. also llj). The slip- link was 
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Figure 4: The loop size probability distribution p{£) ai p — 
55% area coverage, for the F8 with 512 (top line) and 1024 
monomers (solid line). The power law with the predicted 
exponent c = 1.375 in equation Q is indicated by the dashed 
line. 



represented by three tethered beads enforcing the shding 
pair contact such that the loops cannot fuUy retract, see 
figure m To create the dense F8 initial condition, a free 
F8 is squeezed into a quadratic box with hard walls. This 
is achieved by starting oS from the free F8, surrounding 
it by a box, and turning on a force directed towards one 
of the edges. Then, the opposite edge is moved towards 
the centre of the box, and so on. During these steps, the 
slip- link is locked, i.e., the chain cannot slide through 
it, and the two loops are of equal length during the en- 
tire preparation. Finally, when the envisaged density is 
reached, the hard walls are replaced by periodic bound- 
ary conditions, and the slip-link is unlocked. For all den- 
sities studied we observe that one of the loops becomes 
much smaller than the other one. We consider the chain 
as relaxed when the structure has changed from symmet- 
ric (si — S2 — N/2) to asymmetric (si ~ 1, S2 ~ iV), 
then to symmetric and back again to asymmetric. Af- 
ter this we start to sample the probability distribution 
for the loop lengths. The probability distribution with 
chain length N = 1024 and density p — 55% required 
about 3 X 10^ Monte Carlo steps (^^ 3 x 10^^ attempted 
moves). We note that an explicit simulation to obtain 
the PDF p{£) for more complicated structures than the 
F8 was not possible within reasonable computation time 
with the MC algorithm we used. 



IV. THE 2D TREFOIL: SCALING ANALYSIS & 
SIMULATIONS 



We now turn to the 2D projection of the trefoil (de- 
noted "3" , cf. left part and position I in figure^top row). 
Each of the three crossings is replaced with a vertex of 
order four, resulting in a network with Af = 6 seg ments 
of lengths {s;} and total length L = Y^^=i E3- The 



size of the knot region is ^ = J2i=i ^^'^ sq = L — £. Ac- 
cording to section^ we need to know the behaviour of 
the configuration number lu3{£, L) for ^ ^ i in order to 
decide how {£) scales with L. For fixed segment lengths 
{si} with £ L, the configuration number of the net- 
work can be derived by using equation (|A3() in appendix 
^and an analogous reasoning invoked to obtain equation 
resulting in 



co',^coo{L)£''^w[^. 



Si 



£<^L 



(8) 



The prime on lu'^ indicates that the segment lengths {si} 
are fixed. However, since in our case the individual seg- 
ments of the knot region may exchange length with each 
other and only the total length £ of the knot region is 
fixed, we integrate over all distributions of lengths 
{si} under the constraint £ = X]^=i 0- This yields 
the desired number of configurations of the 2D trefoil 
with fixed £ ^ L as 



(9) 



u;3{£,L)^ujo{L)£- 



with c = —73 — TO, where 73 — —33/8 from equation (|A2|) 
in appendix IXI (£ — 4, n4 — 3) and to = 4 is the number 
of independent integrations over chain segments. Thus, 
c = 1/8 < 1 which implies that the dense 2D trefoil is 
delocalised. 

The above analysis corresponds to the case for which 
all segments {si} are large compared with the short- 
distance cutoff a. Conversely, if some of the segments 
{si} are of the order of a, the vertices they join can 
no longer be resolved on macroscopic length scales, but 
constitute a new, single vertex, possibly with more than 
four outgoing legs. The corresponding contractions of 
the original 2D trefoil thus represent different networks, 
each one with its own topological exponent 7 and locali- 
sation exponent c — —7 — to (where m is the number of 
independent integrations, see above). These contractions 
correspond to different shapes, which can be ordered ac- 
cording to their scaling order in a, i.e., according to in- 
creasing powers of c (cf. section [HJ. 

For dilute polymers, we have shown in this way that the 
leading scaling order is the F8 with c = 43/16, and that 
the original (uncontracted) trefoil shape is only found at 
the third position |T^ . For dense polymers, however, 
the present scaling results show that both the original 
trefoil shape (c = 1/8 < 1, see above) and position II 
(c = 3/4 < 1) are in fact delocalised (top row in figure 
The F8 is only found at the third position and is 
weakly localised (c — 11/8 > 1, cf. section ITlT|) . Thus, in 
a MC simulation of the dense 2D trefoil, we predict that 
one mainly observes delocalised shapes corresponding to 
the original trefoil, and less frequently the other shapes 
of the hierarchy (top row) in figure ^ 

These predictions are consistent with the numerical 
simulations of Orlandini et al. [s^l (OSV), who observe 
that the mean value of the second largest segment of the 
simulated 2D dense trefoil configurations grows linearly 
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with L, and conjecture the same behaviour also for the 
other segments, corresponding to the delocahsation of the 
trefoil obtained above. However, we note that the con- 
figuration shown in figure 1 of OSV corresponds in our 
modelling to the shape at position II in figure ^ where 
we consider the two crossings to the right in figure 1 of 
OSV as one "molten" vertex. 

The same reasoning can be applied to the 2D trefoil 
in the 8 phase. In this case we find that the leading 
shape is again the original (uncontracted) trefoil, with 
c — 5/7 < 1. This implies that the 2D trefoil is delo- 
calised also at the Q point. All other shapes are at least 
weakly localised, and subdominant to the leading scaling 
order represented by the original trefoil. The resulting 
hierarchy of shapes is shown in figure ^ (middle row). 

This finding is at variance with the simulation results 
of OSV at the Q point, who observe a behaviour of the 
simulated trefoil configurations similar to the F8, which 
is weakly localised with c = 11/7 > 1, and is found only 
at the third position in the hierarchy of figure ^ 



in 2D one has dv > \ and jo-jvl increases with A'', so 
that the term in the square brackets in the above ex- 
pression increases with A^ as well. For a fixed number 
V = X]Af>4 "Af of vertices this implies that cg is minimal 
if only vertices with four outgoing legs are present. Using 
this and the fact that for A^ = 4 the term in the square 
brackets is negative 40], we conclude that cg is minimal 
if the number of such 4-vertices is maximal. This implies 
that the leading scaling order of the 2D projection of any 
prime knot is the original (uncontracted) configuration, 
for which the above vertices represent the crossings, and 
this configuration will be delocalised (since the 2D trefoil 
configuration is already delocalised). 

We thank J. Cardy and U. Seifert for helpful discus- 
sions. This work was supported by the National Sci- 
ence Foundation (DMR-01-18213) and the US-Israel Bi- 
national Science Foundation (BSF) grant No. 1999-007. 
RM acknowledges partial support from the Deutsche For- 
schungsgemeinschaft . 



V. CONCLUSIONS 

We presented a scaling analysis for the 2D figure-eight 
structure (F8) and the 2D projection of the trefoil knot 
at equilibrium, both in the dense phase and at the O 
point Figure n] shows the hierarchy of contractions 
of the original trefoil, arranged according to scaling or- 
der. The F8 structure represents the leading order of the 
hierarchy in the dilute phase but does not play a 
special role in the other states. Thus we conclude that 
the 2D trefoil is delocalised both in the dense phase and 
at the Q point, in contrast to its localisation in the di- 
lute case [l^ . The delocahsation of the flat trefoil in the 
dense phase has been observed in the simulations by Or- 
landini et al. 'sSj (OSV). However, their observation of 
behaviour corresponding to a weakly localised F8 at the 
Q transition contradicts our results. 

If the chain's topology is that of a F8, then one of the 
loops is predicted to be tight in all cases, although the 
localisation exponent depends on the polymer phase. We 
explicitly verified the localisation exponent of c = 11/8 
by a MC study of a F8 in the dense phase. (We em- 
ployed periodic boundary conditions for this simulations 
to avoid any potential problems associated with surface 
effects.) 

In 2D, the scaling analysis can be readily extended to 
general prime knots. Indeed, the minimal 2D projection 
of any prime knot can be mapped on a network Q, for 
which one can calculate the corresponding exponent c 
in a similar way as before. Using the Euler relations 

SA/" = EAr>i NriN and C = J2n>i li^ ~ '^^n + 1, we 
find 



Cg = 2 + njv 
N>4 



N 



\(Tn\ - dv) 



(10) 



Appendix A: SCALING IN GENERAL POLYMER 
NETWORKS 



In this appendix we review the scaling results for poly- 
mer networks in the dense phase in 2D 25, 26, 27, 2^ 
and at the O transition 35]. 

A general polymer network Q, like the one depicted in 
Fig. [SI consists of a number of vertices which are joined 
by M chain segments of total length L. First, consider 
the dense phase in 2D. If all segments have equal length 
s and L = TVs, the confi gura tion number ug of such a 
network scales as 



he conh gurc 



(Al) 



where ujq{L) is the configuration number of a simple ring 
of length L. For dense polymers, and in contrast to the 
dilute and O phases, ujq{L) (and thus cug) depends on 
the boundary conditions and even on the shape of the 
system |2^ |23 • For example, for periodic boundary con- 
ditions (which we focus on in this study) corresponding 
to a 2D torus, one finds ujo{L) ~ /i-^ with a connec- 

tivity constant fi and 'i' — 1 26]. However, the network 
exponent 



7e 



l-C 



E 

Af>l 



(A2) 



Both for dense polymers and polymers at the O point 



is universal and depends only on the topology of the 
network by the number C of independent loops, and by 
the number njsi of vertices of order N with vertex ex- 
ponents CTAT = (4 - iV2)/32 miH]. For a linear chain, 
the corresponding exponent 7iin = 19/16 has been ver- 
ified by numerical simulations |26l |36j. For a network 
made up of different segment lengths {si} of total length 
L = X^i^i ■^i' equation (|A1|) generalises to (cf. section 4 
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Figure 5: Polymer network Q with vertices (•) of different 
order (m = 5, ns = 4, m = 3, ns = 1). 

in reference [2^) 

ujg{si, ...,sj^) - tJo(i) yg i — , • ■ • , ) , 

\s^^ sj^ J 

(A3) 

which involves the scaling function yg. 

For polymers in an infinite volume and endowed 
with an attractive interaction between neighbouring 
monomers, a different scaling behaviour emerges if the 
system is not below but right at the Q point pjSj . In this 
case the number of configurations of a general network Q 
is given by 

uJg{s,, ...,sm)-h'^ sY' f — , . . . , ^) , (A4) 



with the network exponent 

7g; = 1 - dvC + ^ uno'n ■ (A5) 

Af>l 

Over lined symbols refer to polymers at the O point. In 
d^2, 4/7 and = {2 ^N){2N + 1)/A2 

Finally, we note that equation l)A2|l with ctat = (4 — 
7V^)/32 holds for a 2D dense polymer network which 
never intersects (apart from the vertices). However, it 
has been shown recently that a different, so-called 'Gold- 
stone phase' emerges if the strict non-crossing condition 
is relaxed; i.e., if crossings are allowed albeit, disfavoured 
[28l |. It is argued in reference that in this case the 
scaling dimensions Xk of the so-called A;-leg operators 
all vanish. This implies ctn = in equation l|A2|l (see 
equation (13) in 25j and the Euler relations near equa- 
tion (|10|l above), i.e., the asymptotic behavior is similar 
to a network of chains without self-avoiding constraints 
(ideal polymers). The localisation exponent c is then 
simply given by c = — 7g — m — £ — m — 1, resulting 
in c values of -1, 0, 1, 0, 0, 1, 1 for the 7 contractions 
shown in figure ^ top row from left to right. As ar- 
gued before, all contractions with c < 1 are delocalised; 
in particular the original (uncontracted) trefoil shape, is 
the dominant form with the smallest c. Our finding for 
the 2D dense phase in section Hvl is thus also applicable 
in the Goldstone phase. However, in a MC simulation of 
the Goldstone phase it will be hard to identify the knot, 
as additional crossings are allowed and will occur. 
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